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ABSTRACT 

We calculate the spectrum resulting from the interaction of a fireball with ambient 
soft photons. These photons are assumed to be produced by the walls of a funnel in 
a massive star. By parameterizing the radial dependence of the funnel temperature 
we calculate the deceleration of the fireball self-consistently, taking into account the 
absorption of high energy 7-rays due to interaction with the softer ambient photons. 
The resulting spectrum is peaked at energies in agreement with observations, has a 
v 2 slope in the X-ray band and a steep power-law high energy tail. 

Key words: gamma rays: bursts — X-rays: general — radiation mechanisms: non- 
thermal 



1 INTRODUCTION 

We have recently proposed (Lazzati et al. 2000, hereafter Pa- 
per I) that the gamma-ray burst (GRB) phenomenon orig- 
inates from the interaction of a relativistic fireball with a 
dense photon environment, leading to Compton drag. On 
one hand this is an inevitable effect if the progenitors of 
GRBs are massive stars which are about to explode or 
. have just exploded as supernovae; on the other hand this 
mechanism greatly alleviates the efficiency problem faced 
by the standard internal shock scenario (Lazzati, Ghisellini 
& Celotti 1999; Panaitescu, Spada & Meszaros 1999; Kumar 
1999). In Paper I we have discussed the basic Compton drag 
scenario, showing how this process can convert bulk motion 
energy directly into radiation with a remarkable high effi- 
ciency and, on the basis of simple estimates, how the re- 
sulting spectrum should peak, in a uF(y) representation, 
around ~1 MeV, as observed. 

Here we quantitatively and self-consistently estimate 
the predicted spectrum, assuming that the fireball propa- 
gates in a funnel inside a massive star, and show that, in- 
dependently of the details of the model, it satisfactorily re- 
sembles what observed. Since the funnel walls emit a black- 
body spectrum and the scattered photons are boosted by 
the square of the Lorentz factor (r) of the fireball, the local 
spectrum has a blackbody shape, at a temperature enhanced 
by r 2 . However, the observed spectrum, convolution of all 
the locally emitted spectra, is not a blackbody, due to four 
main effects: i) the funnel walls would not be at a uniform 
temperature, but there should be a gradient between the in- 
ternal and external parts; ii) if the Compton drag process 



is efficient, the fireball decelerates; iii) the very high energy 
emission produced in the internal regions can interact with 
the ambient photons, producing electron-positron pairs; iv) 
the fireball may become optically thin to scattering outside 
the funnel, where the ambient photons are characterized by 
the same temperature, but their energy density is progres- 
sively diluted with distance. 



2 BASIC ASSUMPTIONS 

We postulate that the fireball propagates with a bulk 
Lorentz factor Y inside a funnel cavity, whose walls emit 
blackbody radiation at a temperature T, of conical shape 
with semi-aperture angle ij). The calculation starts at the 
distance zo, assumed to be the end of the acceleration phase 
and, for consistency, we verify that the power emitted at 
2 < zo is negligible. 

We assume that the fireball is and remains cold in the 
comoving frame. At zq, in fact, the internal energy has been 
already used to accelerate the fireball, and thus protons are 
sub-relativistic. On the other hand leptons might be still 
hot at zo and/or being re- heated when the bulk scattering 
process starts to be efficient. However, in a few (Compton) 
cooling timescales they would reach the (sub-relativistic) 
Compton temperature. It is thus reasonable to treat also 
the leptonic component as cold in the estimate of both the 
dynamics and resulting spectrum. 

The initial Lorentz factor (at zo) is indicated as To, 
and the fireball energy is therefore Ef = ToM/c 2 , where 
Mf is its rest mass. For simplicity, the dependence of the 
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temperature on z, between zo and the radius of the star z», 
has been parameterized by a power law: 



T(z) = T 



(i) 



where T„ is the temperature at the top of the funnel. 

Inside it, we approximate the local radiation energy 
density of the ambient photons as U(z) — aT 4 (z). Beyond 
z,, and in the region where the fireball remains optically 
thick (i.e. for z < zt, see below), U(z) is characterized by 
the same temperature, but decreases. As the relevant quan- 
tity is the amount of radiation which is indeed scattered 
by the fireball, we parameterize the dependence on z of the 
product U(z)x (the scattering rate) as (z/z»)~ 9 . 

We consider g a free parameter. Inside the funnel g — 0, 
while outside it a value g > 2 can account for a decrease in 
the scattering rate due to the changing of the typical scat- 
tering angle (photons come preferentially at smaller angles 
as z increases). As the scattering rate is oc (1 — 13 cos 9), 
where 9 is the angle between the photon and the electron 
directions, far from the star surface (1 — /3cos#) oc (z/z»)~ 2 , 
corresponding to g ~ 4. Furthermore some of the radiation 
produced by the massive star could be reflected and re- 
isotropized by scattering material, of unknown radial den- 
sity profile, likely surrounding the massive star progenitor. 
In particular if this forms a wind with a z~ 2 profile, the 
energy density of the re-isotropized radiation scales as z~ 3 , 
and dominates the seed photon distribution at large dis- 
tances. In this case U(z)x (the scattering rate) can have a 
complex profile, being flat in the vicinity of the surface of 
the star, then decreasing as z~ 2 and as z~ 4 for increasing z, 
to become flatter when the component associated with the 
re-isotropized photons dominates. It is also possible that, 
as a result of intermittent stellar activity, the stellar wind 
is not continuous. In this case a single shell may dominate 
the scattering, producing a homogeneous and isotropic scat- 
tered radiation field, dominating the total radiation energy 
density beyond some critical distance. 

The distance zt at which the fireball becomes optically 
thin to scattering is 



zt 



axEf 



nip 2 m p c 2 To 



1/2 



3.7 x lO^-i^ro^'cm (2) 



where the conventional representation Q = Q X W X and c.g.s. 
units are adopted. It is then likely that the fireball becomes 
transparent at z > z* (since the radius of red supergiants is 
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As long as the fireball is opaque to scattering, the inter- 
action with photons boosts their energy by a factor ~ 2T 2 . 
Therefore the (local) total energy emitted by the fireball 
through the Compton drag process (over a distance dz) is 



dE(z) = 2-Kip 2 z 2 aT i (J^j 



dE(z) = 2-K^ 2 z 2 aTt y — 



Y 2 dz 



T 2 dz 



z < z* 



Z > Zt. 



(3) 



(4) 



The factor 2 in front of the RHS of these equations takes 
into account that the preferred scattering angle is ~ 90°, 
corresponding to an average energy boost of 2Y 2 . 

Let us now consider the spectral shape. For this it is 
convenient to use dimensionless photon energies and tem- 



peratures, defined as x = hv/(m e c 2 ) and O = kT/(m e c 2 ), 
respectively. 

The resulting Compton spectrum has a blackbody 
shape, of effective temperature T c = 2r 2 T (or 6 C = 2F 2 9), 
i.e. the local spectral distribution produced within dz is 
given by: 

dE(z,x) = TrV^meC 2 ^) -JL— dz - 

Z < 2* (5) 



s 2 .2Z*(z/z,) a 2 fm e c\ 
dE(z,x) = ■nip p rn e c ^ — j 



-dz: 



e x/e c ,, _ i 
z > z», (6) 

where O c ,» = 2r 2 €>*. Equations (5) and (6) are correctly 
normalized, i.e. the integrated energies correspond to those 
expressed in (3) and (4). 



3 THE FIREBALL DYNAMICS 

As long as the fireball remains optically thick for scattering 
and this occurs in the Thomson regime, the dynamics (de- 
celeration) of the fireball due to the radiative drag, obeys: 

M f c 2 ^- = -2^ 2 z 2 aT i Y 2 . (7) 

Assuming the temperature profile of equation (1) we obtain: 

To 



r = 



l + 27rV> 2 aT ( ?r 2 z3[(z/zo) 3 



■«-l]/[B/(3-46)]' 

zo < z < z», (8) 



and thus the deceleration radius, z<j, defined as the distance 
at which T is halved, corresponds to: 

1/(3-46) 



Zd 



Zo 



1 + 



1 + 



S/(3-4b) 
2nip 2 aT i T 2 l z-> 

E f (3-4b) 



1 1/(3-46) 



2-K^aTiY 2 (z t /z y b 4 



(9) 



Beyond Zd, the Lorentz factor decreases with distance as a 
power law, whose slope is determined by the temperature 
profile. 

Outside the star radius (z > z*) the Lorentz factor fol- 
lows: 



r = 



r, 



1 + 2TnP 2 aT?r r t zt[(z/z t )3-9 - l]/[E f (3 - g)} ' 

z, < z < z T . (10) 



Note that Klein-Nishina effects are important for incoming 
photon energies such that xY > 1, i.e. when > l/(3r). 
For simplicity, we neglect interactions in this regime when 
calculating T(z), but we assume no scattering events when 
Q > l/(3r) in calculating the spectrum. This simplification 
is justified as long as most of the fireball energy is lost in 
the Thomson scattering regime (see Fig. 2, which shows that 
T starts to decrease at distances where the temperature is 
small enough to ensure scatterings entirely in the Thomson 
regime) . 

When the fireball becomes optically thin, the amount of 
scattered photons is correspondingly reduced, and the pro- 
cess becomes less efficient. As shown by equation (2), this 
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is likely to happen at some distance from the star surface, 
where the photon density is also reduced, thus further de- 
creasing the efficiency of the process. In the numerical calcu- 
lations we have however included the optically thin scatter- 
ing regime, and one can see its contribution in Fig. 1 (dotted 
line) . 



4 PAIR PRODUCTION 

A further effect which may strongly affect both the observed 
spectrum and the dynamics of the fireball is the production 
of electron-positron pairs through photon-photon interac- 
tions. Let us thus consider in turn the role of scattered and 
funnel radiation as seed photons for this process. 

4.1 Interaction among photons in the beam 

The threshold energy for interaction between photons of en- 
ergies x and xt is xt > 2/[x(l — cos#)] ~ 4T 2 fx, where all 
quantities are calculated in the observer frame. The latter 
expression takes into account that the high energy photons 
produced by the Compton drag are highly collimated, within 
a typical angle sin 9 ~ 1/F. As the bulk of the scattered 
photons have energies x ~ 2r 2 (30), pair production would 
occur if TO > 1/3. 

However this also implies that the scattering process is 
in the Klein Nishina regime, and we can therefore conclude 
that photon-photon collisions among photons in the beam 
can only affect the high energy tail of the spectrum produced 
at each radius, while the emission at the peak is unaltered. 
We therefore neglect this effect. 




e [keV] 

Figure 1. Example of spectra produced by Compton drag. The 
thick solid lines correspond to the sum of the radiation produced 
inside the funnel (thin solid lines) and outside it (dashed and dot- 
ted lines). The thin solid lines at the highest energies correspond 
to the emission neglecting photon-photon absorption, to show the 
importance of this process. The dotted line (only shown for the 
To = 30 case) is the spectrum produced by the fireball once it 
becomes optically thin. The model parameters are for all cases: 
Ef = 5 X 10 51 erg; i/j = 0.2; b = 0.5; g = 2; z* = 10 13 cm and 
T* = 3 X 10 5 K. The three cases differ for the assumed initial bulk 
Lorentz factor and zo, i.e. To = 30, 100, 300 and zq = 3 X 10 s , 
10 9 , 3 X 10 9 cm, respectively. The two vertical dashed lines mark 
10 and 150 keV, the range of the foreseen hard X-ray detector 
onboard the Swift mission. 



4.2 Interaction between beam photons and funnel 
radiation 

The interaction between the 7-rays produced by the Comp- 
ton drag process and photons emitted by the funnel walls 
would occur at large angles, resulting in an average energy 
threshold xt > 1/x. Since x < To, this absorption mecha- 
nism would be significant as long as the funnel walls produce 
a sufficient number of photons with energies xt > 1/Tq. 

Let us then estimate the photon-photon optical depth 
r 77 , by integrating the product of the photon-photon cross 
section aj 1 (x,xr) and the photon density above threshold 
n-f(x) over the 7-ray path, i.e. from the site of creation, zi, 
to infinity, and over the photon energies: 



dx J <r 77 ( (z, x')dz. 



(11) 



Since cr^(x',x) is peaked at the threshold energy, equation 
(11) can be simplified (Svensson 1984, 1987) as 



f (zi,x) 



5m e c 2 



xtU(z, xr)dz, 



(12) 



where U(z,xt) = m e c?n 1 {z,XT) is the photon energy den- 
sity at threshold, at the location z, i.e. 



U(Z!,X T ) 



8-Tvh 



X* T 



exp[x T /Q{zi)] - 1' 



(13) 



The radiation flux produced at the location z\ is then de- 
creased by the factor exp[— r 77 (zi, x)\ while crossing the fun- 
nel. 

The absorbed radiation will be reprocessed by the pairs, 
and re-distributed in energy. Each electron and positron will 
have an energy 7 ~ x/2 at birth, and will cool due to the 
Compton drag process. The positrons will then annihilate 
in collisions with the electrons in the fireball, producing a 
Doppler blueshifted annihilation line at x ~ F. We have ne- 
glected these reprocessing mechanisms, since, as can be seen 
in Fig. 1, the amount of energy absorbed in 7-7 collisions is 
small, amounting to a few per cent at most. 



5 THE SPECTRUM 

The observed total spectrum can be computed by integrat- 
ing equations (3) and (4) over z, taking into account photon- 
photon absorption. The contribution produced within the 
star is given by: 



E(x) = 



2,2 2(mec\* f z 2 x 3 e-^*> , 
n ^ meC (__) / — dz; 

J Zn 



Zq < z < z+, 



(14) 



while beyond z„ the number of target photons able to inter- 
act with high energy 7-rays to produce pairs is negligible, 
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Figure 2. The profile of the bulk Lorentz factor T corresponding 
to the cases shown in Fig. 1. The vertical dotted line marks 10 13 
cm, the top of the funnel. 
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Figure 3. Spectra produced by Compton drag for three different 
choices of the temperature of the surface of the massive star (as 
labeled). For all cases E f = 5x 10 51 erg; if) = 0.2; b = 0.5; g = 2; 
z, = 10 13 cm, T = 100 and z = 10 9 cm. 



and thus, ignoring photon-photon absorption, we obtain: 

3 r ZT z 2 (z/ Zt )' 3 x 3 



E(x) = jrVW^y 



Y6( e x/Bc,: 

Z* < Z < ZT- 



-dz: 



(15) 



In Fig. 1 we show three examples of the predicted spectrum 
corresponding to different values of the initial bulk Lorentz 
factors. To illustrate the main features of the model and 
the importance of photon-photon absorption, this is calcu- 
lated both with and without the photon-photon absorption 
term. Together with the total spectrum, the separate con- 
tributions for z < 2* and for zt < z < z* are reported. In 
Fig. 2 we show the corresponding T profiles. The effect of 
the star surface temperature (and of the entire funnel, since 
the parameter b is assumed to be the same for all cases) can 
be clearly seen in Fig. 3. Note the v~ x ^ 2 power law shape in 
the X-ray band for the high temperature case. The exten- 
sion of this power law branch depends on the value of g. In 



the case shown (g — 2) the radiation energy density outside 
the funnel remains sufficiently large to cause the decelera- 
tion of the fireball, and this is responsible for the power law 
tail between 10 and 100 keV. For larger g the extension of 
this power law would decrease. This effect can also be seen 
for the high Fo case in Fig. 1. 

In order to determine the general features of the pre- 
dicted spectrum and thus assess its robustness against the 
parameters of the model, we also derived analytical (al- 
though approximated) expressions for the spectral energy 
distribution. 



5.1 Analytical approximations 

First, let us approximate the blackbody spectral form with 
its Ray leigh- Jeans part, and let us neglect photon-photon 
absorption. In this case, for x < 60F 2 we have: 



dE(z,x) oc 



T 



z dz; 



for z < z < z* 



tx pT^ 2 9 dz; for z„ < z < zt- 



(16) 



Three regimes occur at different distances: 
zo < z < zj: — in this case Y = const, and integration over 
z yields: 



E(x) tx x- {3 - 3b)/b ; for z > z d 



(17) 



which, for b = 0.5, gives E(x) tx x 3 . 

Zd < z < z„: — here V decreases as (z / zo) - ' 3-46 ' and thus: 

-3(1-6)7(6-76) iQIZd<z<Zt . 



E(x) oc x 



(18) 

which, for b = 0.5, results in E{x) oc x~ 3/l5 . 
z* < z < zt: — at these distances the ambient radiation 
energy density decreases as (z/zo)~ 9 . If V remains constant 
(= r»), the spectrum E(x) oc x 2 , while, for T decreasing as 

r oc (z/z )" (3 " 9) 



E(x) oc x 



-1/2 



for z t < z < zt , 



(19) 



which is independent of g. 

In conclusion, in the case of efficient Compton drag, 
and independently of the particular choice of parameters, 
the predicted spectrum is always characterized (in order of 
decreasing energy) by: a steep high energy tail; a first break 
flagging the deceleration of the fireball; a second break cor- 
responding to radiation produced at the top of the funnel 
- above which the temperature of the ambient photons re- 
mains constant; a third break, below which the spectrum 
oc x~ x / 2 , corresponding to the deceleration of the fireball due 
to the isothermal photon bath; and finally a fourth break, 
below which the spectrum F(x) oc x 2 . One obtains such a 
hard spectrum, instead of the familiar slope F(x) oc x corre- 
sponding to scatterings of isotropically distributed electrons 
and seed photons, because only the photons scattered along 
the forward direction are observed Q 



* This can be seen by integrating Eq. 7.23 of Rybicki < 
(1979), in the angle range [0 < 6»i < 1/T] 



Lightman 
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6 DISCUSSION 

If the fireball propagates in a dense photon environment the 
Compton drag effect must necessarily be taken into account, 
and it may even be the dominant emission mechanism, able 
to decelerate the fireball without the need of internal shocks 
and without invoking the build-up of large magnetic fields. 

In this letter we have shown that the predicted spec- 
trum, rather than being simply a black body spectrum 
boosted in energy, has a complex shape, with power law 
segments corresponding to the decrease in temperature of 
the funnel, deceleration of the fireball, and dilution of the 
radiation energy density as the fireball propagates outside 
the funnel while remaining optically thick. 

The general features of the predicted spectrum quali- 
tatively agree with observations, since they can explain the 
steep power law high energy tail, the peak of the emission, 
and a hard tail in the X-ray band. The latter feature is par- 
ticularly interesting, since other models made different pre- 
dictions. In the standard internal shock synchrotron model, 
in fact, the spectrum cannot be harder than u 1 ^ 3 in the 
thin part, and it is very unlikely that self-absorption can 
take place in the X-ray band (Granot, Piran & Sari 2000). 
This would in fact imply a huge density of relativistic par- 
ticles, making the inverse Compton effect largely dominate 
the total radiation output. This radiation would be emitted 
at higher and yet unobserved frequencies, and would then 
worsen the already severe efficiency problem. 

In the quasi-thermal Comptonization model, on the 
other hand, the typical predicted spectral shape in the X- 
ray band is oc v° , down to the typical frequencies of the seed 
soft photons, i.e. the IR-optical band (Ghisellini & Celotti 
1999; Meszaros & Rees 2000). 

The existing observations of a significant fraction of 
burst spectra harder than u 1 ^ 3 (Preece et al., 1999a, b; Crider 
et al., 1997) are therefore already a challenge to existing 
models, and may suggest a Compton drag origin of this por- 
tion of the spectrum. However the situation is not already 
a clear-cut because, to receive enough photons to study 
the spectral shape, integration times are much longer than 
the dynamical time-scales of the system, with the spectrum 
rapidly evolving in time. More sensitive instruments, such 
as the Burst Alert Telescope (BAT, a coded mask detec- 
tor more sensitive than BATSE) onboard the foreseen Swift 
satellite will probably overcome this limitation. 

We must also stress that the Compton drag scenario is 
not alternative to the more conventional internal shock one. 
Indeed, the front of the fireball will decelerate first, plausibly 
causing subsequent un-decelerated parts to shock even if the 
central engine is working in a continuous way. This would 
produce additional radiation, either by the synchrotron 
and inverse Compton processes or by quasi-thermal Comp- 
tonization, depending on the details of the particle acceler- 
ation mechanism (see Ghisellini & Celotti 1999). We then 
expect spectral evolution: since the latter radiation mech- 
anisms produce a steeper low energy tail, a hard-to-soft 
transition (i.e. from u 2 to i/ 1 ^ 3 or u°) would occur. 

In this paper, we have considered the illustrative case of 
a single fireball moving out through an extended stellar en- 
velope, along a funnel which is empty of matter but pervaded 
by thermal radiation from the funnel walls. The fireball it- 
self (for typical parameters) remains optically thick until it 



expands beyond the stellar surface. A burst with complex 
time-structure could be modeled by a series of fireballs or 
expanding shells. However, in this more general case, the 
later shells would suffer less drag, since not enough time 
may have elapsed to replenish the entire funnel cavity with 
seed photons. Indeed one expects the spikes to be more pow- 
erful the longer is the time interval between them, as more 
seed photons could pervade the cavity. This, besides causing 
internal shocks with the first shell which has been efficiently 
decelerated by Compton drag, will also result in a distribu- 
tion of T-factors: they will become greater on axis, where 
few seed photons can efficiently Compton drag the shells, 
and smaller towards the border of the funnel, where seed 
photons can be replenished by the funnel walls. 

We plan to investigate these possibilities and their con- 
sequences on the associated predicted afterglows in future 
work. 



ACKNOWLEDGMENTS 

AC acknowledges the Italian MURST for financial support. 



REFERENCES 

Crider A. et al., 1997, ApJ, 479, L39 
Ghisellini G., Celotti A., 1999, ApJ, 511, L93 
Granot J., Piran T., Sari R., 2000, ApJ, 534, L163 
Kumar P., 1999, ApJ, 523, L113 

Lazzati D., Ghisellini G., Celotti A. & Rees M.J., 2000, ApJ, 529, 
L17 (Paper I) 

Lazzati D., Ghisellini G. & Celotti A., 1999. MNRAS, 309, L13 
Meszaros P. & Rees M.J., 2000, ApJ in press, |astro-ph/9908126 
Panaitescu A., Spada M. & Meszaros P., 1999, ApJ, 522, L105 
Preece R.D. et al., 1998a, ApJ, 496, 849 
Preece R.D. et al., 1998b, ApJ, 506, L23 

Rybicki G.B. & Lightman A. P., 1979, Radiative Processes in As- 
trophysics, J. Wiley & Sons (New York) 
Svensson R., 1984, MNRAS, 209, 175 
Svensson R., 1987, MNRAS, 227, 403 



© 0000 RAS, MNRAS 000, 000-000 



